7 5AJ iooob 


NASA TECHNICAL NOTE 


■ 

o 



<x, 

oo 


EXTRACTION FROM FLIGHT DATA OF 
LATERAL AERODYNAMIC COEFFICIENTS FOR 
F-8 AIRCRAFT WITH SUPERCRITICAL WING 

by James L. Williams and William T. 

Langley Research Center 
Hampton, Va. 23665 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • NOVEMBER 1974 






1. Report No. 2. Government Accession No. 

NASA TN D-7749 

3. Recipient's Catalog No. 

4. Title and Subtitle 

extraction from flight data of lateral 

AERODYNAMIC COEFFICIENTS FOR F-8 AIRCRAFT 
WITH SUPERCRITICAL WING 

5. Report Date 

November 1974 

6. Performing Organization Code 

7. Author(s) 

James L. Williams and William T. Suit 

8. Performing Organization Report No. 

L-9619 

10. Work Unit No. 

501-26-04-02 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


16. Abstract 


A parameter -extraction algorithm was used to determine the lateral aerodynamic 
derivatives from flight data for the F-8 aircraft with supercritical wing. The flight data 
used were the recorded responses to aileron or rudder pulses for Mach numbers of 0.80, 
0.90, and 0.98. 


Results of this study showed that a set of derivatives were determined which 
yielded a calculated aircraft response almost identical with the response measured in 
flight. Derivatives extracted from motion resulting from rudder inputs were somewhat 
different from those resulting from aileron inputs. It was found that the derivatives 
obtained from the rudder -input data were highly correlated in some instances. Those 
from the aileron input had very low correlations and appeared to be the more reliable. 


17. Key Words (Suggested by Author(s)) 

18. Distribution Statement 

Parameter extraction 

Unclassified - Unlimited 

Maximum likelihood 

Lateral aerodynamic coefficients 

STAR Category 02 


19. Security Oassif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22, Price* 

Unclassified 

Unclassified 

35 

$3.25 


For sale by the National Technical Information Service, Springfield, Virginia 22151 





















EXTRACTION FROM FLIGHT DATA OF 


LATERAL AERODYNAMIC COEFFICIENTS FOR F-8 AIRCRAFT 
WITH SUPERCRITICAL WING 

By James L. Williams and William T. Suit 
Langley Research Center 

SUMMARY 

A parameter -extraction algorithm was used to determine the lateral aerodynamic 
derivatives from flight data for the F-8 aircraft with supercritical wing. The flight data 
used were the recorded responses to aileron or rudder pulses for Mach numbers of 0.80, 
0.90, and 0.98. 

Results of this study showed that a set of derivatives were determined which yielded 
a calculated aircraft response almost identical with the response measured in flight. 
Derivatives extracted from motion resulting from rudder inputs were somewhat different 
from those resulting from aileron inputs. It was found that the derivatives obtained from 
the rudder -input data were highly correlated in some instances. Those from the aileron 
input had very low correlations and appeared to be the more reliable. 

INTRODUCTION 

Analytical and simulator studies of the flight and handling qualities of aircraft 
require that accurate estimates of the aerodynamic parameters be used if the results are 
to be valid. To provide aerodynamic parameters for analytical and simulator studies and 
also to provide numerical values for comparison with wind-tunnel data and theoretical 
estimates, parameters have been extracted from flight data for many years. Results 
from recent Langley Research Center investigations in which the aerodynamic parameters 
for several aircraft have been extracted from flight data are reported in references 1 to 5. 

Various techniques for data extraction are presented in reference 6. The technique 
and program used in extracting the parameters in this study is the maximum likelihood 
method of reference 7. 

The purpose of the present study is to extract static and dynamic lateral aerody- 
namic parameters from flight tests of the F-8 supercritical -wing aircraft at Mach 
numbers of 0.80, 0.90, and 0.98. The flight-test data used for the analysis consisted of 
airplane responses to aileron or rudder pulses. 


SYMBOLS 


a Y 

acceleration measured along Y body axis, g units 

b 

wing span, m 

DA = 5 a - 

«a,t» rad 

DR = 6 r - 

6 r ,t, rad 

g 

gravitational acceleration, 9.81 m/sec^ 

h P 

pressure altitude, m 

Jx^Y^Z 

moment of inertia about roll, pitch, and yaw axis, respectively, kg-m 

*XZ 

product of inertia, kg-m^ 

*h 

distance from aircraft center of gravity to quarter -chord point of 
mean geometric chord of horizontal tail, m 

M 

Mach number 

N 

number of data points 

P 

roll rate, rad/sec 

q 

pitch rate, rad/sec 


R 1 

a Y W ay 

R -J- 

Rp 'w p 


W v 


2 


2 


yaw rate, rad /sec 



S 


u 

V 

v 

W 

W aY ,W p ,W 
W V ,W y 


wing area, 

velocity component along 
resultant velocity, m/sec 
velocity component along 
weight, N 

weighting values of 



X body axis, m/sec 


Y body axis, m/sec 


a y , P, r, v, and 


(p 


state variables 


w velocity component along Z body axis, m/sec 

X,Y,Z body-axis system (see fig. 1) 

x vector of aircraft states 

a angle of attack, deg or rad 

/3 angle of sideslip, deg or rad 

6 a aileron deflection, total deflection of left and right ailerons, positive 

deflection indicates positive roll, deg or rad 

6 r rudder deflection, positive trailing edge left, deg or rad 

6 pitch attitude, deg or rad 

p mass density of air, kg/m 3 

(p roll attitude, rad 
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Subscripts: 
c computed 

m measured in flight 

t trim 

Aerodynamic parameters: 

rolling-moment coefficient about 




ac, 

Qjg, = -gf Per radian 

dCj 

Ci. = — - per radian 

96 a 


= -zir P er radian 


c n yawing -moment coefficient about 



ac n 

C n/3 = P er radian 


9C n 

C n = per radian 
n 5 a 8 5a 


body axis 


body axis 
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C n 5r . = W Per 

Cy lateral -force coefficient along Y body axis 



9Cy 

C v = — - per radian 

Y 0 9/3 

3Cy 

C v = per radian 

Y 6 a 96 a 


9Cy 

C-u- = — — per radian 

Y 6 r 36 r 


The dot over a symbol denotes rate of change with respect to time. 


EQUATIONS OF MOTION 


The equations of motion used in this study are referred to a body -axis system 
(fig. 1) and are as follows: 

Y-direction: 
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Yawing: 


r = -pq 


h~ l X 

iz 


far - » T? + 5 + V s ' *> + C "P H 


+ ^ nr 2V + ^ n 6 a ^ a " + ^ n 5 r ^ r " §r >^ 


Auxiliary equations: 
a = tan- 1 J 

fi = Bin’ 1 £ 


1 /• 

ay = — (v + ur - wp - g cos 6 sin <p) 


( 3 ) 


The three-degree-of -freedom lateral equations (eqs. (1) to (3)) were solved for the 
lateral motions. The longitudinal state variables u, w, and q used in these equations 
were the flight -measured values, and hence the nonlinear contributions of these terms 
were used in the equations. 


DESCRIPTION OF AIRPLANE 

The aerodynamic parameters were extracted from data obtained from an F-8 air- 
craft on which the original wing was replaced with a supercritical wing which is described 
in reference 8. The F-8 is a single -seat high performance airplane with a single jet 
engine embedded in the fuselage and with an all-movable horizontal tail. Figure 2 is a 
photograph of the airplane. Pertinent geometric characteristics of the airplane are 
given in table I. 

Flight instrumentation appropriate to this study included the following items: 

(1) Roll-, pitch-, and yaw -rate gyros 

(2) Angle -of -attack indicator 

(3) Angle -of -sideslip indicator 

(4) Altimeter 

(5) Total velocity indicator 

(6) Aileron and rudder position indicators 

(7) Magnetic tape recorders 
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(8) Accelerometers 

(9) Pitch- and roll -angle gyros 

The full-scale range of the flight instruments and their accuracies are given in 
table II. 


FLIGHT TESTS 

The data presented in this study were obtained from in-flight measurement of the 
airplane response to aileron or rudder inputs. These flights were made at the NASA 
Flight Research Center as part of a general evaluation program of the F-8 supercritical- 
wing aerodynamics. 

Data were obtained at Mach numbers of 0.80, 0.90, and 0.98. The control inputs 
used to generate the lateral motions (three degrees of freedom) were one or more aileron 
or rudder pulses. All data used in this study for parameter extraction were reduced at 
Flight Research Center and have been corrected for bias and displacement of the measur- 
ing instrument from the aircraft center of gravity. Test conditions and mass character- 
istics are presented in table III. 

PARAMETER-EXTRACTION PROCEDURE 


The parameter -extraction procedure used in this study is an iterative technique 
which utilizes the maximum likelihood method to extract the stability and control param- 
eters. (See ref. 7 for details of the method.) This method maximizes the conditional 
likelihood L which is a function of aerodynamic parameters, weights, and initial 
conditions: 


L = 


1 

(2it) 1 /2r1/2 



~ x i,c)^R l( x i,m 



where R is the estimate of the error covariance matrix and x is a vector describing 
the state of the aircraft. Maximization of the likelihood function L provides the follow- 
ing information: 

(1) The parameter increments which are used to update the parameter 

(2) The covariance matrices whose elements are proportional to the estimated 

standard deviations and the pairwise correlation coefficients for the param - 
eters and for the states 
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(3) The performance index function J which is an indicator of the fit between 
measured and calculated motion 

The interative technique produces a set of estimated derivatives which, when used 
in the equations of motion (eqs. (1) to (3)), provide the best fit to the time variation of the 
aircraft motion measured in flight. The criterion for the best fit is the performance 
index J which is defined as 


J - 



x i,c)( x i,m 


* 1 ,c) T 


where det denotes determinant and x is the vector describing the state of the air- 
craft. Generally, the performance index J becomes smaller with sucessive iterations. 
The iteration procedure is stopped when the value of J does not change appreciably for 
several successive iterations. 

The components of the vector x are the state variables v, p, r, cp, and a Y 
which are calculated from equations (1) to (3). The quantities u m , v m , and w m were 
not measured directly but were obtained from the measured total velocity, angle of side- 
slip, and angle of attack through the use of the following equations: 


u m = v m cos «m cos &n 
v m - V m sin 0 m 
w m = V m cos j3 sin a m 

Maximization of the likelihood function yields the covariance matrix for the meas- 
urement noise based on the current solution. (See ref. 7.) This matrix gives the vari- 
ances (or standard deviations) of the differences of the measured state and the current 
solution. The inverse of this matrix is the weighting matrix used in the parameter - 
change equations. In this investigation, the diagonal form of the weighting matrix was 
used and the diagonal elements can be expressed as the squares of the differences 
between the measured and calculated data. For example, the weight for the state vari- 
able v is expressed as 


I (vm-vc)f 

R V 1=1 
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Similar equations are obtained for the other state variables used in the cost function 

(y^ 'Ml 'Mir™* vH) 

Initial values of the state variables were obtained from the flight records for the 
time period prior to a control input. Initial values of the aerodynamic derivatives 
(table IV) were obtained from a preliminary investigation in the Langley 8 -foot transonic 
pressure tunnel for a = 4°. Although the tabulated derivatives were for a Mach number 
of 0.8, they were used as initial values for all the work presented herein. The data rate 
for the in-flight measurements was 25 points per second. 

RESULTS AND DISCUSSION 

Measured and calculated time histories are compared in figures 3 to 6. The cal- 
culated motions are those obtained after the parameter -extraction program had converged 
on a set of derivatives for each flight condition. Figures 3 and 4 show the responses to 
rudder inputs for flight at a Mach number of 0.80 and 0.90, respectively. Figures 5 and 6 
are responses to aileron inputs for a Mach number of 0.90 and 0.98, respectively. In all 
cases, the calculated time histories match the flight results very closely. The standard 
deviations of the converged states are given in table V. 

The aerodynamic parameters extracted for each of the four flight conditions in 
figures 3 to 6 are listed in table VI. Also shown in the table are the estimates of the 
standard deviations of the parameters. The standard deviations are generally less than 
10 percent of the extracted derivative, which implies that the extracted derivatives are 
well defined. Most of the derivatives are of reasonable magnitude and are close to the 
initial values shown in table IV. However, a comparison of derivatives for the rudder - 
input and the aileron-input cases at the same Mach number (M = 0.90) shows appreciable 
differences in some derivatives, for example, C^, Cj , and C^. The difference in 
C; r was particularly bothersome since there was a sign difference. In order to deter- 
mine if this discrepancy was caused by differences in control input or for some other 
reason, the correlation matrices for each set of extracted derivatives were examined. 

The correlation matrices are shown as tables VII (aileron input) and VIH (rudder input). 

It is apparent that several high correlations exist fo** the rudder -input case, whereas the 
correlations are quite low for the aileron -input case. It appeared, therefore, that the 
derivatives from the aileron-input time histories were more reliable than those from the 
rudder -input cases. In order to check this possibility, the aerodynamic derivatives that 
had been extracted from the rudder -input flight data were used to calculate motions for 
the aileron-input case. The results are presented in figure 7 and show poor agreement 
with flight data. Motions were then calculated for the rudder -input case from the aero- 
dynamic parameters extracted from the aileron-input data. The results, shown in 
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figure 8, show good agreement with the flight data. The conclusion, therefore, is that 
the derivatives extracted from the aileron -input data are more reliable than those from 
the rudder-input data. The precise reason for the superiority of an aileron input over 
a rudder input in obtaining good, uncorrelated aerodynamic parameters is not known at 
this time, and is a subject for future study. As a matter of interest, the derivatives 
from the aileron -input case were used as initial values for the rudder -input flight (which 
resulted in the comparison of fig. 8), and the parameter -extraction program was per- 
mitted to iterate. After several iterations the fit to the data improved, and the extracted 
parameters converged to the values of the derivatives of table VI for the rudder input. 

The negative values of Cj r extracted from the aileron-input data, theoretically, 
should have been positive. It was believed that perhaps the time histories were insensi- 
tive to Cjr r ; therefore, several computations were made for the rudder -input case with 
M = 0.9 for all derivatives fixed at the values of table VI except Cj , The results are 
shown in figure 9, and indicate that the effect of varying Cj r is quite noticeable. 

Some of the extracted parameters of table VI are compared with wind-tunnel values 
in table IX. In general, the agreement is reasonably good; however, there are some 
appreciable differences between the flight -extracted and wind-tunnel values. The most 
obvious difference is the large variance between the extracted and wind-tunnel values for 
C n6 . One reason for this is possibly aeroelastic effects. Since rigid body equations of 
motion were used in the extraction algorithm, the effect of aeroelasticity on the extracted 
parameters was not considered in this study. 

CONCLUDING REMARKS 

A parameter -extraction algorithm was used to determine the lateral aerodynamic 
derivatives from flight data for the F-8 aircraft with supercritical wing. The flight data 
were the responses to aileron or rudder pulses for Mach numbers of 0.80, 0.90, and 0.98. 

Results of this study showed that a set of derivatives were determined which yielded 
a calculated aircraft response almost identical with the response measured in flight. 
Derivatives extracted from motion resulting from rudder inputs were somewhat different 
from those resulting from aileron inputs . It was found that the derivatives obtained from 
the rudder -input data were highly correlated in some instances. Those from the aileron 
input had very low correlations and appeared to be the more reliable in fitting the meas- 
ured responses. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., August 12, 1974. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF F-8 AIRCRAFT 
WITH SUPERCRITICAL WING 


Fuselage: 

Length, m 16.09 

Wing: 

Area, m 2 25.45 

Aspect ratio 6.77 

Span, m 13.14 

Mean geometric chord, m . . 2.08 

Vertical tail: 

Area, m 2 10.13 

Aspect ratio 1.5 

Span, m > . . . 3.88 

Rudder: 

Area, m 2 1.17 

Horizontal tail: 

Area, m 2 8.68 

Aspect ratio 3.5 

Span, m 5.52 

Tail length, center of gravity to quarter-chord point of 

mean geometric chord, m . . . .' 5.31 
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TABLE n.- INSTRUMENT RANGES AND ACCURACIES 


Measured quantities 

Range 

Accuracy 

Angle of attack, deg '. 

-5 to 25 

±0.3 

Pressure altitude, m 

0 to 18 000 

±11 

Total velocity, m/sec 

0 to 360 

±0.7 

Lateral acceleration at c.g., g units 

±1.0 

±0.01 

Longitudinal acceleration at c.g., g units 

±0.5 

±0.01 

Roll angle, deg 

±90 

±1.8 

Pitch angle , deg 

±30 

±0.6 

Roll rate, deg/sec 

±200 

±0.6 

Pitch rate, deg/sec 

±40 

±0.8 

Yaw rate, deg/sec 

±12 

±0.3 

Rudder position, deg 

±8 

±0.3 

Aileron position, deg 

±15 

±0.8 


TABLE HI.- TEST CONDITIONS AND MASS CHARACTERISTICS 



M = 0.80 

M = 0.90 

M = 0.98 

Weight, N 

100 802 

104 912 

101 792 

Mass, kg 

10 279.1 

10 698.2 

10 380.0 

I X , kg-m 2 

20 484 

20 512 

20 491 

Iy, kg-m 2 

124 383 

125 350 

124 616 

Iz, kg-m 2 

138 424 

139 363 

138 650 

l XZ> k g-m 2 ' 

4359 

4522 

4399 

h p , m ' . . . 

12 178 

12 751 

13 862 

*h- m 

5.31 

5.30 

5.27 
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TABLE IV.- INITIAL VALUES OF AERODYNAMIC DERIVATIVES 


% 

% 

C, 





-1.317 

-0.234 

-0.390 

0.069 

-0.492 

-0.005 

0.231 

-0.183 

0.148 


TABLE V.- STANDARD DEVIATIONS OF CONVERGED STATES 


■ M 

Input 

Standard deviation of - 

v, 

m/sec 

P, 

rad/sec 

r, 

rad/ sec 

<P, 

rad 

a Y , 

g 

0.80 

Rudder 

0.4095 

0.0036 


0.0018 

0.0036 

.90 

Rudder 

.3041 

.0065 


.0045 

.0046 

.90 

Aileron 

.3080 

.0085 


.0091 

.0044 

.98 

Aileron 

.3326 

.0165 

.. 

.0028 

.0163 

.0067 
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TABLE VI.- EXTRACTED AERODYNAMIC PARAMETERS AND STANDARD DEVIATIONS 


Parameters 

Parameters and standard deviations a for - 

Rudder input 

Aileron input 

M = 0.80 

M = 0.90 

M = 0.90 

M = 0.98 

C Y,t 

-0.01597 

-0.0172 

0.0014 

-0.0171 


(0.000149) 

(0.00013) 

(0.00017) 

(0.00017) 

c y r 

-1.2063 

-1.2055 

-1.2283 

-1.1868 

P 

(0.009819) 

(0.0086) 

(0.0209) 

(0.02363) 

C Y 6 

0.03182 

0.0325 

0.0320 

0.0320 


(0.001093) 

(0.00127) 

(Fixed) 

(Fixed) 

C Y 6a 

0 

0 

-0.0228 

-0.01342 


0 

0 

(0.00357) 

(0.00325) 

c i;t 

-0.000202 

0.00051 

-0.00053 

-0.00034 


(0.000012) 

(0.000013) 

(0.00002) 

(0.00002) 

% 

-0.23399 

-0.2052 

-0.2748 

-0.2751 


(0.00200) 

(0.00174) 

(0.0018) 

(0.00284) 

% 

-0.4948 

-0.3928 

-0.5938 

-0.6867 


(0.00614) 

(0.00529) 

(0.0044) 

(0.00743) 

c *r 

0.2111 

1.5767 

-0.4546 

-0.4749 


(0.03553) 

(0.03256) 

(0.0342) 

(0.0887) 

% 

0.0046 

0.0063 

0.005 

0.005 

o r 

(0.000069) 

(0.000075) 

(Fixed) 

(Fixed) 

C 4 a 

0.0940 

0.0940 

0.0941 

0.0740 


(Fixed) 

(Fixed) 

(0.00051) 

(0.00062) 

^■n,t 

-0.000093 

-0.00028 

0.00022 

0.00000 


(0.000007) 

(0.000009) 

(0.00001) 

(0.000009) 

C n^ 

0.1449 

0.1552 

0.1473 

0.1334 


(0.000901) 

(0.00091) 

(0.00041) 

(0.00069) 

C„ 

0.02292 

0.00796 

-0.0059 

-0.0360 

y 

(0.002674) 

(0.0026) 

(0.00166) 

(0.00311) 

C n r 

-0.20581 

-0.31961 

-0.4368 

-0.6057 


(0.02029) 

(0.0230) 

(0.01051) 

(0.0241) 

C n g 

-0.01481 

-0.01526 

-0.015 

-0.0150 

°r 

(0.000047) 

(0.000059) 

(Fixed) 

(Fixed) 

Cn 6 a : * ’ * 

-0.0022 

-0.0022 

-0.0022 

-0.0024 


(Fixed) 

(Fixed) 

(0.00019) 

(0.00026) 


a Standard deviations are given in parentheses. 
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TABLE VII.- ESTIMATED PAIRWISE CORRELATION COEFFICIENTS AT CONVERGENCE FOR AILERON INPUT AT M = 0.90 












































































































































TABLE VHL- ESTIMATED PAIRWISE CORRELATION COEFFICIENTS AT CONVERGENCE FOR RUDDER INPUT AT M = 0.90 



































































































































TABLE IX.- COMPARISON OF EXTRACTED DERIVATIVES WITH 
WIND-TUNNEL RESULTS 


Coefficient 

M = 

0.80 

M = 0.90 

M = 

0.98 

Wind 

tunnel 

Rudder 

input 

Wind 

tunnel 

Rudder 

input 

Aileron 

input 

Wind 

tunnel 

Aileron 

input 

C, .... 
L V 

-0.390 

-0.495 

-0.420 

-0.393 

-0.594 

-0.445 

-0.687 

c „ r — 

-.492 

-.206 

-.583 

-.320 

-.437 

-.601 

-.606 

% ' ‘ ‘ ' 

.148 

.145 

.166 

.155 

.147 

.206 

.133 

c Y/3 ... 

-1.327 

-1.206 

-1.450 

-1.206 

-1.228 

-1.501 

-1.187 

C h ' ‘ ‘ * 

-.234 

-.234 

-.304 

-.205 

-.275 

-.172 

-.275 

C« .... 
n P 


.023 


.008 

-.006 

-.010 

-.036 

c \ • • ■ 

.069 


.069 


.094 

.069 

.074 

c n 6r • • • 

-.183 

-.015 

-.183 

-.015 


-.183 


























Time, t, sec 


Figure 3.- Comparison of flight data with time histories computed by using the 
aerodynamic parameters of table VI for rudder input at M = 0.80. 
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Figure 5.- Comparison of flight data with time histories computed by using the 
aerodynamic parameters of table VI for aileron inputs at M = 0.90. 
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Figure 7.- Concluded. 
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